Momentum Dependence of Resonant Inelastic X-Ray Scattering Spectrum in 

Insulating Cuprates 
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The resonant inelastic x-ray scattering spectrum in insulating cuprates is examined by using the 
exact diagonalization technique on small clusters in the two-dimensional Hubbard model with second 
and third neighbor hopping terms. When the incident photon energy is tuned near the Cu K 
absorption edges, we find that the features of the unoccupied upper Hubbard band can be extracted 
from the spectrum through an anisotropic momentum dependence. They provide an opportunity 
for the understanding of the different behavior of hole- and electron-doped superconductors. 



PACS numbers: 74.25. Jb, 71.10.Fd, 78.70.Ck 

Resonant inelastic x-ray scattering (RIXS) is develop- 
ing very rapidly into a powerful technique to investigate 
elementary excitations in the strongly correlated electron 
systems [^j-Q . The application of this technique to insu- 
lating copper oxides has made it possible to observe an 
excitation due to a local charge transfer between cop- 
per and oxygen and local d-d excitations on copper 
site [||. In addition, it has been demonstrated that, by 
using high resolution experiments the momentum- 
dependent measurement of the charge transfer gap is 
possible when the incident photon energy uji is tuned 
through Cu K absorption edge. Thus, the RIXS can be 
a useful probe to obtain information on the momentum 
dependence of the elementary excitations. 

One of the elementary excitations in the insulating 
cuprates is the charge-transfer process from the occu- 
pied Zhang-Rice singlet band (ZRB) |^ composed of 
Cu 3dx^-y^ and O 2p^ orbitals to the unoccupied up- 
per Hubbard band (UHB). The dispersion of ZRB have 
been extensively studied by angle-resolved photoemission 
spectroscopy (ARPES) experiments on the parent com- 
pounds of high Tc superconductors |^-|^ : A d- wave-like 
dispersion was observed along the (0,7r)-(7r,0) line with 
the minimum of the binding energy at (7r/2,7r/2) On 
the contrary, the dispersion relation and spectral proper- 
ties of the unoccupied UHB have not been examined and 
thus remain to be understood. The information of UHB 
is of crucial importance for the understanding of the mo- 
tion of electrons in the electron-doped superconductors. 
In addition, it may be useful to know if the particle-hole 
symmetry is required for the high temperature supercon- 
ductivity. 

In this Letter, we examine the RIXS spectrum for the 
Cu K-edge, and demonstrate that the characteristic fea- 
tures of the dispersion of UHB can be extracted from 
the momentum dependence of the spectrum. To see this, 
we use the half-filled single-band Hubbard model to de- 
scribe the occupied ZRB and unoccupied UHB by map- 
ping ZRB onto the lower Hubbard band (LHB) in the 
model. Then, we incorporate Cu Is and 4p orbitals into 
the model to include the Is-core hole and excited 4p elec- 



tron into the intermediate state of the RIXS process. The 
long-range hoppings are also introduced in the Hubbard 
model with realistic values obtained from the analysis 
of ARPES data. We find a characteristic momentum 
dependence of the Cu K-edge RIXS spectrum: The en- 
ergy of the threshold of the RIXS spectrum at (7r/2,7r/2) 
is higher than that at (0,0), whereas the energy of the 
threshold at (7r/2,0) is lower than that at (0,0). This 
anisotropic dependence is explained by the dispersion of 
the UHB which has the minimum energy at (7r,0) due 
to the long-range hoppings. The determination of the 
UHB will contribute to the understanding of the differ- 
ent behavior of hole- and electron-doped superconduc- 
tors 1,0. 

We map the ZRB onto the LHB, which is equivalent 
to the elimination of O 2p orbitals. Such mapping was 
used in the analysis of O Is x-ray absorption spectrum 
pd| . The Hubbard Hamiltonian with second and third 
neighbor hoppings for the Sd electron system is written 
as, 

-i" E 4X- + H.c. + f/E<T"U' (1) 

(ij)3rd,0- i 



where d\^ is the creation operator of 3d electron with 

spin a at site i, nf^ — dj^dj^rj, the summations (i,j)ist, 
(i7j)2nd, and (i,j)3rd run over first, second, and third 
nearest-neighbor pairs, respectively, and the rest of the 
notation is standard. 

Figure ^ shows the schematic process of Cu i^T-edge 
RIXS. An absorption of an incident photon with energy 
LUi, momentum K^, and polarization e,; brings about the 
dipole transition of an electron from Cu Is to 4p orbital 
[process (a) in Fig. ||. In the intermediate states, 3d elec- 
trons interact with a Is-core hole and a photo-excited Ap 
electron via the Coulomb interactions so that the excita- 
tions in the 3d electron system are evolved [process (b)]. 
The 4p electron goes back to the Is orbital again and a 
photon with energy wy, momentum Kj, and polarization 
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FIG. 1. Schematic picture of the Cu K-edge RIXS process. 
An incident photon is absorbed and dipole transition Is 4p 
is brought about [process (a)], and through the intermediate 
state [process (b)], the photo-excited ^giicWsost^& iib Is 
again and a photon is emitted [process (c)]. 

£/ is emitted [process (c)] . The differences of the energies 
and the momenta between incident and emitted photons 
are transferred to the 3d electron system. 

In the intermediate state, there are a Is-core hole and 
a 4p electron, with which 3d electrons interact. Since the 
Is-core hole is localized because of a small radius of the 
Cu Is orbital, the attractive interaction between the Is- 
core hole and 3d electrons is very strong. The interaction 
is written as, 



H 



ls-3d 



= -V 



"■i,cr"i,(T' 



(2) 



where is the number operator of Is-core hole with 
spin G at site i, and V is taken to be positive. On the 
contrary, since the 4p electron is delocalized, the repul- 
sive interaction between the 4p and 2>d electrons as well 
as the attractive one between the 4j) electron and the Is- 
core hole is small as compared with the Is-Sd interaction. 
In addition, when the core-hole is screened by the id 
electrons through the strong Is-Sd interaction, effective 
charge that acts on the 4p electron at the core-hole site 
becomes small. Therefore, the interactions related to the 
4p electron are neglected for simplicity. Furthermore, we 
assume that the photo-excited 4p electron enters into the 
bottom of the 4^^ band with momentum ko, where z-axis 
is perpendicular to the Cu02 plane. This assumption is 
justified as long as the Coulomb interactions associated 
with the 4p electron are neglected and the resonance con- 
dition is set to the threshold of the Xs^^pz absorption 
spectrum [Q. Under these assumptions, the RIXS spec- 
trum is expressed as. 



/(AK,Aa;)=^ 



Sko-K/,o-Pko,o- 
1 



,-p?'ko,o-'*ko-Ki,o-l'^) 



H — Eq — UJi — iT 

xd{Auj - Eo, + Eq), 



where H = Had + i?is-3d + -ffis,4p, ^^15,4^ being kinetic 
and on-site energy terms for a Is-core hole and a 4p elec- 
tron, AK = Ki - K/, Aw = - ujf, s^^^ (Pk,cr) is the 
creation operator of the Is-core hole (4p electron) with 
momentum k and spin a, |0) is the ground state of the 
half- filled system with energy Eq, \a) is the final state of 
the RIXS process with energy Ea , and F is the inverse of 
the relaxation time in the intermediate state. In Eq. (|^), 
the terms His^^p are replaced by eis_4p which is the en- 
ergy difference between the Is level and the bottom of 
the Apz band. 

The RIXS spectrum of Eq. (||) is calculated on 
(VSxVS)-, {VTOxVW)-, and (4x4)-site clusters with 
periodic boundary condition by using a modified ver- 
sion of the conjugate-gradient method together with the 
Lanczos technique. We will show the results for the 4x4- 
site cluster in the following. 

The values of the parameters are as follows: t'/t = 
-0.34, t"/t = 0.23, U/t = 10, V/t = 15, and T/t = 1 
with t = 0.35 cV. The values of t, t' , and t" are the same 
as those used in the analysis of ARPES data of the high 
Tc superconductors based on the t-t'-t"-J model [p|,p^. 
The value of U is obtained from the relation J=4i^/[/ 
and J/t—OA. The value of V is set to be larger than 
that of U. The results shown below are insensitive to the 
magnitude of V as well as of F |l^ . 

Before going into the RIXS spectrum, we mention the 
resonance condition in RIXS. To determine the condi- 
tion, we have to examine the Cu Is x-ray absorption 
spectroscopy (XAS) spectrum defined as. 



D{uj) = -Im(0|sko-K,,aPko,^77 ^ ^ 



XPko,a4o-K„.|0), 



(4) 



(3) 



where H is the same as that in Eq. (||) . It is necessary to 
tune the incident photon energy uji to the energy region 
where the Cu Is XAS spectrum appears. The inset in 
Fig. H shows -D(w), where a two-peak structure appears, 
i.e., the broad one around Lu—Sis-4:p = —20t and the sharp 
one around —13t. The former mainly contains configura- 
tions that the core-hole site is doubly occupied by the 3c? 
electrons {U — 2V = — 20i), while the latter dominantly 
contains configurations that the core-hole site is singly 
occupied {~V — — 15t). This means that, when the in- 
cident energy is tuned around the former structure, the 
information about UHB can be extracted from the RIXS 
spectrum. Thus, we set oji to the threshold of the XAS 
spectrum denoted by the arrow in the inset. 

Figure]^ shows the momentum dependence of the RIXS 
spectrum. The spectra below Auj/t ~ 2 and above 
Au/t ~ 5 have different origins: The former comes from 
the excitations related to the spin degree of freedom such 
as two-magnon Raman scattering, the energy scale of 
which is so small that the spectrum is hard to be ob- 
served M ■ On the other hand, the latter is related to the 
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FIG. 2. Resonant inelastic x-ray scattering spectra for 
Cu -ft'-edge in half-filled Hubbard model with long-range 
hoppings. The spectra of the elastic scattering process at 
AK=(0,0) are not shown. The parameters used are U /t = 10, 
V/t = 15, T/t = 1, t'/t = -0.34, and t" /t ^Fl^-.^TsTfeS ff^^i" 
cal dotted hne denotes the position of the peak at AK — (0, 0) 
for guide to eyes. The J-functions (the vertical thin solid lines) 
are convoluted with a Lorentzian broadening of Q.2t. Inset is 
the Cu Is absorption spectrum with FxAs/t = F/f = 1.0, and 
the incident photon energy uji is set to the value denoted by 
the arrow. 



excitations from LHB to UHB. The vertical dotted line in 
the figure denotes the position of the low-energy peak at 
AK = (0, 0) for guide to eyes. The spectra strongly de- 
pend on the momentum showing a feature that the weight 
shifts to higher energy region with increasing |AK|. This 
momentum dependence is also obtained in the 10-site 
cluster calculations. In addition, the threshold of the 
spectrum at AK = (7r/2,0) and (7r,7r/2) is lower in en- 
ergy than that at (0,0). At AK = (7r/2,7r/2), how- 
ever, the spectrum appears above the threshold at (0,0), 
resulting in an anisotropic momentum dependence be- 
tween the spectra along (0,0) to (7r/2,0) and along (0,0) 
to (77/2, 7r/2). We note that such an anisotropic feature 
can not be observed in the RIXS spectrum of the Hub- 
bard model without t' and t" (not shown) In the 
8-site cluster the spectrum at AK = (7r/2,7r/2) is lower 
in energy than that at (7r,0), which is consistent with the 
above feature. 

In the intermediate state, the excitations of the 3d elec- 
trons from the occupied to unoccupied states are caused 
by the interaction with the Is-core hole, i.e., i?is-3d = 

The operator d\^^ (j'^k2,<T represents a particle- hole exci- 
tation. Therefore, we analyze the RIXS process by de- 
composing into such a particle-hole excitation in order 
to understand the anisotropic momentum dependence of 
the threshold along (0,0) to (7r/2,0) and along (0,0) to 
(^/2,7r/2) in Fig. |. 
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FIG. 3. Single particle excitation spectrum A(]i,uj) in the 
half-filled Hubbard model with t' and t" terms. t'Jt = ^0.34 

Fig. 3, Isiitsui el al. 

and t" jt — 0.23. The dotted line denotes the chemical po- 
tential. The 5 functions are convoluted with a Lorentzian 
broadening of 0.2t. 



As a first step, we consider the particle-hole excita- 
tion as the convolution of the single-particle excitation 
spectra A(k, uj) between occupied LHB and unoccupied 
UHB. Figure g shows A(k, w) in the half-filled Hubbard 
model with t' and t" terms Below the chemical 

potential denoted by the dotted line, a sharp peak ap- 
pears at (7r/2,7r/2) with the lowest-binding energy. In 
contrast, the spectrum at (7r,0) is very broad and deep in 
energy. These features are consistent with the ARPES 
data for Sr2Cu02Cl2 [@,||. Above the chemical potential, 
the dispersion of UHB has the minimum of the energy at 
k=(7r,0) |lj]. Wc show below that the (7r,0) spectrum 
in UHB plays a crucial role in the RIXS spectrum with 
AK=(7r/2,0). 

Now, we examine the lowest-energy excitations with 
AK=(0,0), (7r/2,0), and (7r/2,7r/2) in the convoluted 
spectrum J^^^dEACk + AK,£ + Ljj)A{'k,£), where /i is 
the chemical potential. For the case that AK=(0,0), the 
minimum excitation energy in the convoluted spectrum 
is ^5t at k=(7r,0). In the same way, the lowest-energy 
excitation with AK=(7r/2,0) in the convoluted spectrum 
is that from k=(7r/2,0) of LHB to (7r,0) of UHB with 
the energy of ~4<. This value is smaller than that for 
the AK=(0,0) case, being consistent with the relation of 
the thresholds of the RIXS spectra between AK=(0,0) 
and (7r/2,0). In contrast, the lowest-energy excitation 
with AK=(7r/2,7r/2) in the convoluted spectrum is in- 
consistent with that in the RIXS spectrum, because the 
excitation energy, which is determined by the peaks at 
k=(7r/2,-7r/2) of LHB and at (7r,0) of UHB, is almost 
the same as that for AK=(7r/2,0). This means that the 
argument based on the convoluted spectrum of A(k, lu) 
is insufficient to understand the anisotropic behavior of 
the threshold, and thus we need to treat the process of 
the particle-hole excitation exactly. Therefore, we intro- 
duce the spectral function i?(k, AK;u;) of the two-body 
Green's function, which can describe the particle-hole ex- 
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FIG. 4. Spectral function i3(k,AK,aj) which repre- 
sents the particle-hole excitation process f f 16*©!"' <■' j(anel) , 
with [k,AK] = [(7r/2,0), (7r/2,0)] (solid line) and 
[(7r/2, — 7r/2), (7r/2, 7r/2)] (dotted line). Upper panel shows 
the corresponding RIXS spectra /(AK, Aoj) which are given 
in Fig. |. 



citation, defined as 



i3(k,AK;c.) = ^ 



iE4 



k+AK,(T"k,(T| 



(5) 



where the states \a) have the same point-group symmetry 
as that of the final states of the RIXS process. Figure ^ 
shows S(k,AK;w) with [k, AK] = [(7r/2, 0), (7r/2, 0)] 
and [(7r/2, -7r/2),(7r/2,7r/2)]. Note that k -I- AK = 
(tt, 0) in both cases. The spectrum with [k, AK] = 
[(77/2, 0), (7r/2, 0)] reproduces very well the RIXS spec- 
trum near the threshold. The intensity of _B(k, AK;a;) 
with AK — (7r/2,7r/2) is very small, implying that the 
process X]<t ^(,r,o),o-^('^/2 -'r/2),o- is almost forbidden [|8). 
This is the origin of the small weight of /(AK,a;) with 
AK = (7r/2, 7r/2) around uu/t ^ 5. We note that this be- 
havior can not be obtained by the convolution of A(k, uj) 
mentioned above. At the higher energy region, ui/t ^ 7, 
different processes keeping AK — (7r/2,7r/2), for exam- 
ple, the annihilation of the (7r/2, 0) electron and the cre- 
ation of the (tt, 7r/2) one, dominate the excitations. These 
processes induce the large weight above uo/t ^ 7. 

In summary, we have examined the momentum depen- 
dence of the Cu K-edge RIXS spectrum by using numer- 
ically exact diagonalization technique on small clusters. 
Regarding the ZRB as the LHB, we have adopted the 
Hubbard model with Cu Is- and 4p-bands. We have also 
introduced the long-range hoppings, t' and t'\ of the ZR 
singlet, and found that the threshold of the spectrum 
at AK = (7r/2,0) is small compared with that at (0,0), 
whereas the threshold at (7r/2,7r/2) is larger than that 
at (0,0). This dependence is caused by the (7r,0) state in 
UHB. Thus, by examining the RIXS spectrum, we can 
extract the property of the unoccupied states that is cru- 
cially important for the electron-doped superconductors 
and also for the understanding of the behavior differ- 
ent from the hole-doped superconductors [p[pX|]. Very 



recently, Stanford's group has reported [|I9[ that inter- 
esting data with momentum dependent inelastic scatter- 
ing signal has been seen in the 2 eV range, which re- 
veals the property of unoccupied states right above the 
charge-transfer gap. These progress in both theory and 
experiment will open a new prospect of the physics of 
cuprates. 
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